Noncommutative QCD corrections to the gluonic decays of heavy quarkonia 
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We compute the Noncommutative QCD (NCQCD) contributions to the three gluon decay modes 
of heavy quarkonia. For triplet quarkonia (ortho-quarkonia), the NCQCD correction to the QCD 
three gluon decay mode, like the standard model contribution, is infrared finite. In the case of singlet 
quarkonia (para-quarkonia) , whose QCD three gluon decay mode has infrared singularities which 
are removed using one-loop corrections to the two gluon mode, we find that NCQCD contribution 
is also infrared finite. The calculations are performed in the weak binding limit and do not require 
the introduction of additional effective couplings. 



1. INTRODUCTION 

Efforts to explore the physical implications of field 
theories formulated on noncommutative spaces have 
increased recently 0, Q due to developments in string 
theories ^ |^ , which suggest that noncommutative 
field theories are well defined quantum theories 4] . In 
noncommutative geometry, the coordinates obey 
the commutation relations 

[x^' ,x•']^ie^"' , (1) 

where O'^'^ = —0^^. A noncommutative version of an 
ordinary field theory can be obtained by replacing all 
ordinary products with Moyal -k products defined by 

U^9){x) = e^v{\e^^''^,.^y.)!{x)g{y)l^^ . (2) 

Here, we use the generalization of the QCD Lagrangian 

£ = ^*^V- W*V-7r9Tr(^^p.*F^") , (3) 
where 

D^ii = 9^?/' - * V' , (4) 

and 

\A^, Ai,'\^ = A^ -k Ai, — An -k . (6) 

When supplemented with a gauge fixing term, includ- 
ing a ghost contribution, the Lagrangian Eq. (jS)) can 
be used to obtain a set of vertices and Feynman rules 
for perturbative calculations JS, 7] . 
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2. NCQCD CORRECTIONS TO THE 
LIFETIMES OF HEAVY ORTHO AND PARA 
QUARKONIA 

In standard model QCD, the hadronic contributions 
to the widths of ground state quarkonia are attributed 
to gluonic decays. The pseudoscalar states, para- 
quarkonia, decay predominantly into two gluons while 
the vector states, ortho-quarkonia, being color singlet 
spin one states, must decay into three gluons. Unlike 
para-positronium, which cannot decay into three pho- 
tons in ordinary QED due to charge conjugation sym- 
metry, para-quarkonium can decay into three gluons. 
However, the three gluon mode is infrared singular and 
one-loop corrections to the two gluon mode must be 
included to obtain a finite contribution to the hadronic 
width 1^. 

By assuming that quarkonia are weakly bound it 
is possible to calculate the NCQCD correction to the 
three gluon lifetimes by computing the annihilation 
amplitudes for a noninteracting quark and antiquark 
at rest and supplying a factor of the square of the 
bound state wave function at the origin, 1-0(0)1^, to 
account for the leading binding effect. There is no 
need to devise an effective interaction as in the case of 
the pion decay into three photons (9|,ilfl|. The NCQCD 
amplitudes contributing to the three gluon corrections 
were calculated using the Feynman rules of Ref . 0] • 
These rules contain contributions involving O^n of the 
form k'^6^,yk2, where ki and ^2 are the momenta of 
two gluons. To ensure that the unitarity conditions 
ei^^et^" > and ef^nXpOt^^e^P = are satisfied, we take 
^ofe = —Oka = and write 9ij as 

% = -1^2— 1^13 kOk, (7) 

where 9k is a unit vector and Ajvc is the noncommuta- 
tivity scale. We then have ki9fj,uk2 = 9-{ki x A:2)/A^^. 

The amplitudes to be calculated are illustrated in 
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Fig-ffl Diagrams with the three gluon final state con- 
nected to the quark line by a single virtual gluon do 
not contribute. We computed the threshold helicity 



FIG. 1: The left diagram is illustrative of the class of con- 
tributions to the three gluon decay of ortho-quarkonium 
and the right diagram is one of the class of contributions 
to the para-quarkonium three gluon decay. 

amplitudes with the aid of the symbolic manipula- 
tion program FORM and simplified the results 
using Mathematica 12j. In what follows, we consider 
quarkonium in the "non-relativistic" approximation, 
taking the quark (and the anti-quark) 4-momentum 
p'^ to be p = (m, 0,0, 0). The gluon 4-momenta are 
labelled ki, with i = 1,2, 3. 

All the squared amplitudes can be expressed as func- 
tions of the parameters s,t,u given by 



s = {ki + k2f, ^ = (fc2 + /c3)^ 

which satisfy the identity 

s + t + u ~ 4m^ 



u = (fc3 + fci)' (8) 



(9) 



In terms of these variables the squared matrix element 
summed over the gluon helicities is, for the triplet state 
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— /(5 + 4sin2(^)(s + t + M) (10) 

{t + uf + t^{u + sf + u^{s + tf 
{s + tf{t + uf {u + sf ' 
and for the singlet state 

\^M\^ = —g^{%-Asm^Lp){st + tu + usf (11) 
9 

s^ + t'^ + u'^ + {s + t + uf 

1 ' 

stu{s + tf [t + uf {u + sf ' 



where 



(12) 



The integration over phase space for either state is 
straightforward and the expression for the width in 
terms of the scaled variables x = s/Am?, y = t/Am^ 
takes the form 



dcosS 192(27r)3 



1 l-x 

dx J dy\Mqq^3g\ 




(13) 



The variable tp can be expressed in terms oi x, y and 
the dimensionless scale parameter z = m^/A^^ as 

1, 



Lp = -j/ci xA:2||6'| cos(5 — ^y xy{l — x ~ y)zcosS . (14) 



2.1. Ortho-quarkonium 



The ortho-quarkonium decay width can be sepa- 
rated into two terms, ^^g^^g and F^'^l^^^.The first, 
which is independent of z, gives, after completing the 
phase space integration, the standard QCD result [Tsf 



-.QCD 
oQ^3g 



40 
81 



\m\' 



(15) 



Using the symmetry of the integrand with respect 
to the variables s, t and u (or, equivalently, with re- 
spect to the 4- vectors ki, k2 and ^3) the (z dependent) 
NCQCD contribution can be written 



rpNCQCD 
oQ^3g 

dcosS 



16 ,\mi 

27 ' m2 



1 pI-x 

dx / dy 
Jo 



(16) 



x'^ sin^ ^^/xyy{T^~x~^y) z cos 
(x + y)2(l-y)2 

Integration over dcosS (— 1 < cos (5 < 1) gives 
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NCQCD _ 16 3 1-0(0)1" 



oQ^3g 



27 



.(17) 



2z^/xy{T^^xr^ 



2z^xy(\ -x-y) 



Due to the presence of the square root in the argument 
of the sine function in Eq. (|17|) it is not possible to 
perform the integration analytically. However, rather 
than simply keeping the leading term in z, we inves- 
tigated the behavior of the correction to all orders in 
z. The result, when combined with Eq. (|15|l . has the 
form 



ioQ^3g - -^a. 



5(7r2 - 9) 



+ -(385 - 39^2)^2^(^) 



(18) 



where the behavior of /(z) is illustrated in Fig.|2] 



2.2. Para-quarkonium 

The evaluation of the QCD contribution to Eq. (|13|l 
in the case of para-quarkonium is complicated by the 
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FIG. 2: The function f{z) is appearing in Eq. I18II is plot- 
ted. 



FIG. 3: The function g{z) is appearing in Eq. 1221 is plot- 
ted. 



existence of infrared divergences. By requiring the 
variables x and y to satisfy x,y > e and using the 
symmetry of the integrand, the QCD contribution can 
be written 



^QCD _ 
pQ-*3g ~ 



|V(0)| 



l-2e 

dx 



1 — e—x 

dy 



(19) 



[l + Zx^){x^ -{x + y){l-y)y 



xy{l - x- y){l - xY{l - yY{x + yY 
which can be integrated to give 

(20) 



pQCD _ 3 



2 [152-Il7r2 
41og(e)(ll + 61og(£))] 



The infrared behavior of r^Q_:!^3 exhibited in Eq. ^ 
must be combined with the one-loop corrections to the 
two gluon decay to obtain a finite correction to the 
hadronic width @. 

The evaluation of the non-commutative contribution 
involves an additional factor of —4 sin^ Lp, which may 
be handled as in the ortho-quarkonium case to obtain 



^NCQCD 
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l-2e 

dx 



l~e—x 

dy 



{l + 3x^){x^~{x + y)(l-y)y 
xy{l — X — y){l — xY{l — yY(x -\- yY 

sin [2z^Jxy{\ — x — y) 



2z^Jxy{l -x-y) 



With the introduction of sin^ iy9, the integrand of 
Eq. H21|l is no longer singular when xy{\ — x — y) 
and we may complete the integration with e = 0. The 
NCQCD contribution then can be written 



r 



NCQCD 4 3\jMfi)r_fo-7^2 

81 TO^ 



pQ^Sg 

where g{z) is shown in Fig.|3| 
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z^giz), (22) 



3. DISCUSSION AND CONCLUSIONS 

The inclusion of NCQCD corrections to the three 
gluon decay widths of ortho- and para-quarkonium 
does not change the magnitudes of their hadronic 
widths substantially for scales A^c of order 1 Tev. 
That said, the NCQCD results are interesting in the 
sense of what they imply about the consistency of per- 
turbative calculations in these models. 

For ortho-quarkonium, the two gluon decay is for- 
bidden by Yang's theorem. Hence the three gluon de- 
cay cannot have any soft gluon singularities because 
there are no two gluon one-loop corrections available 
to cancel them. The absence of soft gluon singulari- 
ties, which is, of course, a feature of the QCD three 
gluon decay width, persists when the small (positive) 
NCQCD correction is included. This is guaranteed by 
the form of the squared amplitude, Eq. (|10|l . 

The situation with para-quarkonium is some- 
what more involved because, unlike para-positronium, 
which can only decay into two photons, both two gluon 
and three gluon decays are allowed in QCD. In this 
case, the three gluon decay has infrared singularities, 
which can be combined with the one-loop corrections 
to the two gluon decay to obtain a finite contribution 
to the hadronic width |Si,|l4l- Here, too, the NCQCD 
correction is infrared finite, but only as a result of a 
cancellation provided by the NCQCD effective cou- 
pling. Interestingly, this correction is negative, and, 
(21) while small for realistic values of z = m^/A^^, it can- 
not change the sign of the total width for any value of 
z. 

In summary, we are led to the conclusion that it is 
not necessary to invoke the smallness of z = m? / K\(-, 
to obtain a sensible NCQCD correction to the hadronic 
decays of quarkonia. In principle, should we be pre- 
sented with a fourth generation of very heavy quarks 
with lifetimes sufficiently long to produce quarko- 
nia, corrections to their decay widths associated with 
non-commutative geometry could be calculated con- 
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sistently in the sense that they are finite for any value 
of z. 
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